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Abstract.  The  cyclic  prefix  system  is  used  in  discrete  mul- 
titone  channels  for  frequency  domain  equalization.  This  is 
based  on  the  inversion  of  samples  of  the  channel  frequency 
response  at  uniformly  spaced  points  (DFT  coefficients).  In 
this  paper  we  consider  nonuniformly  spaced  samples  of  the 
channel  frequency  response,  especially  ocatve-spaced.  The 
anticipated  advantage  is  that  for  channels  with  rapidly 
decaying  frequency  responses,  there  are  relatively  fewer 
equalizer  coefficients  with  unduly  large  values,  and  this 
helps  to  reduce  amplification  of  channel  noise  at  the  re¬ 
ceiver.  We  show  how  to  combine  a  wavelet-like  filter  bank 
with  traditional  DFTs  to  achieve  this  goal.  While  the  idea 
appears  to  be  exciting  and  opens  up  interesting  problems, 
its  merits  still  remain  to  be  evaluated.* 

I.  INTRODUCTION 

Discrete  multitone  modulation,  which  today  is  used  in  DSL 
technology,  employs  the  cyclic  prefix  system  for  frequency 
domain  equalization  [7].  The  introduction  of  cyclic-prefix 
with  length  equal  to  the  order  of  the  (FIR)  channel  C(z) 
ensures  that  the  channel  can  be  equalized  by  using  the  in¬ 
verse  of  C(e^u)  at  uniformly  sampled  set  of  frequencies 
—  2nk/M  (Fig.  1(a))  where  M  is  the  block  length 
used  in  the  cyclic  prefixing  [6].  It  is  well  known  that  when 
the  samples  C(e3UJk)  get  very  small,  the  frequency  domain 
equalizers  1  /C{e^k)  amplify  channel  noise  severely.  If  the 
samples  can  be  crowded  more  effectively  near  the  regions 
where  \C(e^)\  is  large  as  in  Fig.  1(b),  this  problem  can 
be  reduced  significantly.  In  this  paper  we  develop  struc¬ 
tures  which  allow  us  to  obtain  such  nonuniform  samples 
by  using  a  combination  of  traditional  uniform  DMT  sys¬ 
tems  and  nonuniform  filter  banks.  While  the  idea  is  quite 
interesting,  further  work  is  required  to  examine  its  practi¬ 
cal  merits.  There  are  a  few  open  problems  that  still  need 
to  be  addressed. 

Preliminaries .  Standard  multirate  notations  (decima- 
tors,  expanders,  delay  chains,  and  so  forth)  will  be  used 
throughout,  and  can  be  found,  for  example,  in  [8].  Thus 
|  M  and  |  M  represent  the  M -fold  decimator  and  ex¬ 
pander  respectively,  and  [X(z)]im  represents  the  z  trans¬ 
form  of  the  decimated  version  x{Mn).  The  blocked  version 
a(n)  of  a  signal  a(n)  with  block  length  K  is  defined  as 
shown  in  Fig.  2,  and  Fig.  3  shows  the  unblocking  opera¬ 
tion.  It  will  be  assumed  throughout  that  the  channel  is  an 


FIR  filter 
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C(z)  =  '£c(n)z~n 
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followed  by  additive  noise  e(n)  (Fig.  4).  The  FIR  assump¬ 
tion  is  usually  justified  because  of  the  channel  shortening 
(or  time-domain  equalization)  step  that  normally  precedes 
cyclic  prefixing  [4]. 
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Fig.  1.  Sampling  the  channel  frequency  response  for  equal¬ 
ization.  (a)  Uniform  sampling,  and  (b)  nonuniform  sam¬ 
pling. 
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Fig.  2.  (a)  The  blocking  schematic,  and  (b)  an  example. 
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Fig.  3.  (a)  The  unblocking  schematic,  and  (b)  an  example. 


The  fact  that  C  is  circulant  implies  that  it  can  be  diago¬ 
nalized  by  the  DFT  matrix  [5].  More  precisely 


noise  e(n) 


c  =  W-1AeW 

where  W  is  the  M  x  M  DFT  matrix  with  elements 
[W]*m  =  Wkm 


Fig.  4.  A  linear  channel  with  additive  noise. 


with  W  =  e~2lri/M.  Here 


II.  THE  TRADITIONAL  CYCLIC  PREFIX  SYSTEM 

Fig.  5(a)  shows  how  a  vector  signal  s(n)  is  cyclic-prefixed. 
Here  s  (n)  is  an  M -component  vector 

So  (*t)  ' 

«i  (ft) 

LsAf-i(n). 


and  p(n)  is  a  copy  of  the  last  L  components  of  s(n),  that 
is, 


p(n)  = 


SM-i(n) 


SM-l(ft) 


The  ( M  +  L)-ve ctor 


Ac  =  diag{C[0],  C[l], . . . ,  C[M  —  1]} 
where  C[k ]  are  the  DFT  coefficients  of  the  channel,  i.e., 
M—l  L 

C[k]  =  ^2  c(n)Wkn  =  Y2c{n)e~i2vkn/M 

n=0  n=0 

Evidently  C[k]  are  samples  of  the  channel  resposne  C{e3<jj) 
at  the  uniformly  spaced  DFT  frequencies  UJk  =  2tt k/M, 
If  these  are  nonzero  for  all  M  values  of  k ,  then  C  is  non¬ 
singular,  and  we  can  equalize  the  effect  of  the  channel.  In 
practice  the  equalization  is  done  by  inserting  the  DFT  and 
IDFT  matrices  as  shown  in  Fig.  6.  The  signal  component 
s(n)  sees  the  cascade  of  W^jC,  and  W  which  is  equal 
to  Ac.  This  is  cancelled  by  the  equalizer  A~*.  If  the  chan¬ 
nel  frequency  response  C(eJu;)  is  close  to  zero  near  any  of 
the  DFT  frequencies,  then  C[k]  is  very  small  there.  This 
means  that  A”1  has  large  numbers  in  it,  and  there  could 
be  severe  noise  amplification. 


x(n)= 


p(n) 

s(n) 


is  said  to  be  the  (L,  M)  cyclic  prefixed  version  of  s (n).  This 
is  then  unblocked  to  obtain  the  scalar  signal  x(n),  which 
is  transmitted  on  the  channel.  The  received  signal  y(fl)  is 
blocked  with  block  length  L  +  M ,  and  the  first  L  compo¬ 
nents  of  each  block  deleted.  The  resulting  M-vector  y(n) 
is  the  received  version  of  the  transmitted  vector  s(n).  The 
usefulness  of  cyclic  prefix  systems  arises  when  the  channel 
can  be  assumed  to  be  FIR  with  order  L  <  M.  In  this  case 
the  transformation  from  s(n)  to  y (n)  is  the  simple  system 
shown  shcematically  in  Fig.  5(b),  where  CisaMxM  cir¬ 
culant  matrix,  and  e(n)  is  the  M-blocked  version  of  the 
channel  noise  e(n).  For  example,  when  M  =  6  and  L  =  3 


e(n) 


Fig.  5.  (a)  The  basic  idea  of  the  (L,  M)  cyclic  prefix 
system,  and  (b)  equivalent  schematic. 


Fig.  6.  The  cyclic  prefix  transceiver. 


The  cyclic  prefix  system  just  described  will  be  denoted  as 
the  (L1M)  cyclic  prefix,  to  remind  ourself  that  the  two 
integers  L  (prefix  length)  and  M  (original  block  length) 
are  involved  in  the  definition.  In  practical  applications  the 
vector  signal  s(n)  is  typically  derived  from  a  scalar  symbol 
stream  in  one  of  two  possible  ways. 


Since  there  are  M  4-  L  samples  in  the  same  duration  that 
was  originally  occupied  by  M  samples,  there  is  a  band¬ 
width  expansion  ratio  7  =  ( M  4-  L)/M  which  can  be 
regarded  as  the  price  paid  for  the  cyclic  prefixing.  This 
means  that  the  symbol  spacing  at  the  channel  input  is 
closer  than  the  original  spacing,  as  schematized  in  Fig. 
8.  Extensions  of  the  cyclic  prefix  technique  in  fractionally 
spaced  equalization,  and  in  blind  equalization  methods  are 
described  in  [9]  and  [10]. 

III.  NONUNIFORM  DMT  SYSTEMS 

We  explain  the  idea  of  nonuniform  DMT  with  an  example, 
which  is  readily  generalized.  The  first  step  in  the  devel¬ 
opment  is  to  split  the  original  binary  symbol  stream  into 
substreams  with  different  symbol  spacings. 


Fig.  7.  The  (B,  M)  parser,  (a)  Parsing  a  binary  sequence 
s(n)  into  M  groups,  B  bits  at  a  time,  and  (b)  schematic 
block  diagram. 


Fig.  9.  The  nonuniform  delay  chain  system,  (a)  Two 
channel  perfect  reconstruction  system,  (b)  analysis  tree, 
and  (c)  synthesis  tree. 


Fig.  8.  Schematic  of  bandwidth  expansion. 

1.  Parsing.  This  is  the  most  common  situation.  Here  a 
scalar  binary  symbol  stream  s(n )  is  divided  into  12- 
bit  blocks,  and  each  block  divided  into  M  groups  as 
demonstrated  in  Fig.  7(a)  for  M  =  3.  The  fcth  group 
then  defines  a  &/-- bit  symbol  stream  Sfc(ra)  (possibly 
converted  to  a  PAM  or  QAM  format)  which  forms  the 
kth  component  of  s (n),  thereby  defining  s(n)  com¬ 
pletely  (Fig.  7(b)).  We  shall  refer  to  the  parser  as 
the  (B,AT)  parser  ( B  bits  turned  into  M  symbol 
streams).  The  inverse  system  at  the  receiver  which 
turns  the  M  symbol  streams  Sfc(n)  back  into  the  bi¬ 
nary  stream  s(n)  is  called  the  (B,  M )  unparser. 

2.  Blocking .  In  some  applications  s(n)  is  a  vector  ob¬ 
tained  by  blocking  a  scalar  symbol  stream  s(n),  pos¬ 
sibly  a  PAM  or  QAM  stream  (Fig.  2). 


Fig.  10.  Equivalent  nonuniform  analysis/synthesis  chain. 


To  accomplish  this,  consider  Fig.  9(a)  which  shows  a  two- 
channel  delay  chain  filter  bank  [8]  with  the  perfect  recon¬ 
struction  property  y(n)  =  x(n).  Suppose  the  two-channel 
analysis  banks  are  connected  in  a  tree  to  obtain  the  system 


of  Fig.  9(b),  and  similarly,  synthesis  banks  connected  as 
in  Fig.  9(c).  This  analysis/synthesis  system  is  equivalent 
to  the  nonuniform  perfect  reconstruction  system  shown  in 
Fig.  10.  In  this  system,  the  input  stream  x(n)  is  merely 
partitioned  into  substreams  of  different  rates,  as  demon¬ 
strated  in  the  figure.  The  synthesis  system  simply  inter¬ 
leaves  them  appropriately  to  obtain  x(n)  again. 


binary 


Fig.  11.  The  nonuniform  DMT  filter  bank,  (a)  Transmitter 
side,  and  (b)  receiver  side. 


Fig.  12.  The  subDMT  systems  in  the  nonuniform  DMT 
filter  bank,  (a)  Transmitter  side  and  (b)  receiver  side. 


Now  consider  Fig.  11(a).  This  is  again  an  example  with 
four  channels,  and  demonstrates  the  main  features  of  the 
system.  The  reciprocals  of  the  decimation  ratios  (2, 4, 8, 8) 
add  up  to  unity.  The  input  binary  stream  is  first  par¬ 
titioned  nonuniformly,  and  then  the  substreams  are  pro¬ 
cessed  by  subDMT  blocks.  Each  subDMT  block  performs 
the  operation  described  in  Fig.  12(a)  for  the  fcth  channel. 


Thus,  there  is  (Bk,Mk)  parsing  of  the  signal  Xjt(n),  fol¬ 
lowed  by  IDFT  and  then  (Lfc,Mfc)  cyclic  prefixing.  The 
resulting  vector  is  then  unblocked  to  obtain  Sjt(n).  The 
subDMT  blocks  create  bandwidth  expansion  by  a  factor 

7k  =•  {Lk  +  Mk)/Mk 

Even  though  the  quantities  Bk,Mk,  and  L *  need  not  be 
the  same  for  all  the  subDMT  blocks,  we  will  assume  that 
7 k  is  the  same  for  all  fc.  This  ensures  that  the  inputs  to 
the  transmitter  filters  F^{z)  in  Fig.  11(a)  have  identical 
symbol  spacing  for  all  k .  The  filter  outputs  can  therefore 
be  added  to  get  the  signal  s(n)  which  is  eventually  trans¬ 
mitted  on  the  channel. 

The  receiver  for  the  nonuniform  DMT  system  is  shown 
in  Fig.  11(b).  Notice  that  the  set  of  filters  (F)t(z),  Hk{z)} 
acts  as  a  (nonuniform)  transmultiplexer  for  the  channel  as 
shown  separately  in  Fig.  13.  The  channel  output  r(n)  is 
first  separated  by  the  receiver  filters  Hk(z)  into  subbands 
rfc(n).  Each  of  the  subband  signals  then  goes  through  the 
receiver  counterpart  of  the  subDMT  system,  as  described 
in  Fig.  12(b):  the  cyclic  prefix  is  first  removed,  then  a 
DFT  performed,  and  finally  the  unparser  generates  the 
received  version  t/fc(n)  of  the  symbol  stream  Xfc(n).  Ev¬ 
idently  yjfc(n)  has  different  spacing  for  different  fc.  Inter¬ 
leaving  these  using  the  nonuniform  delay  chain,  we  finally 
obtain  y(n)  which  is  the  received  version  of  x(n).  Needless 
to  say,  if  the  input  x(n)  in  Fig.  11(a)  is  not  binary  but  an 
arbitrary  signal,  we  can  replace  the  parser  and  unparser  in 
Fig.  12  by  blocking  and  unblocking  operations. 
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Fig.  13.  The  nonuniform  transmultiplexer. 

IV.  OPERATION  WITH  IDEAL  FILTERS 

Continuing  our  example  of  the  four  channel  nonuniform 
system  in  Fig.  11,  let  us  assume  that  the  filters  F ^(z) 
are  ideal  with  responses  as  in  Fig.  14(a).  Assuming  Hk(z) 
have  identical  responses,  it  follows  that  i^fc(eJW)i?m(eJW)  = 
0  for  A;  ^  m  so  that  there  is  no  crosstalk  between  sub¬ 
channels  of  the  transmultiplexer.  The  system  of  Fig.  13 
therefore  merely  splits  the  channel  into  four  subchannels 
with  unequal  bandwidths.  The  path  from  s*(n)  to  r*(n) 
is  the  filter 

Ck{z)  =  [Hk(eju)Fk(^)C(^)]^  (1) 

where  j  njt  is  the  decimator  in  the  fcth  branch.  For  exam¬ 
ple,  the  quantity  Hi(e^)Fi(e^u;)C(e^u;)  is  the  piece  of  the 
channel  highlighted  in  Fig.  14(b)  and  the  8-fold  decimation 


yields  Fig.  14(c).  Thus  the  ideal  nonuniform  transmulti¬ 
plexer  merely  separates  different  portions  of  the  channel 
and  spreads  them  out  into  the  region  [0, 27r].  The  subDMT 
systems  perform  equalization  by  sampling  the  subchannel 
frequency  responses  as  shown  by  the  vertical  lines  in  Fig. 
14(c).  Thusrthe  set  of  samples  of  the  channel  created  by 
all  the  different  paths  appears  nonuniform  as  demonstrated 
in  Fig.  14(d)  which  assumes  DMT  block  length  M  =  4 
for  simplicity  (in  practice  this  could  be  different  from  the 
number  of  transmultiplexer  channels). 
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Fig.  14.  (a)  Example  of  ideal  filter  responses,  (b)  portion 
of  channel  C(eju})  in  the  band  of  Fi(z),  (c)  decimation  by 
8,  and  (d)  nonuniform  sampling  of  the  channel  effected  by 
the  nonuniform  DMT  system. 

We  can  realize  the  nonuniform  trans multiplexer  by  build¬ 
ing  a  tree  structured  filter  bank  with  two-channel  mod¬ 
ules.  Assume  for  example  that  a  two-channel  filter  bank 
costs  N  multipliers  per  symbol.  Then  because  of  the  dec¬ 
imation  involved,  the  total  cost  of  the  tree  structure  is 
N  +  0.5N  +  0.25N . . .  «  2 N.  With  filters  of  moderate  or¬ 
der,  e.g.,  N  =  51,  we  can  get  excellent  approximations  of 
the  ideal  response,  and  cross  talks  created  by  the  transmul¬ 
tiplexer  can  be  ignored.  In  practice  the  subDMT  system 
can  also  choose  not  to  allocate  bits  at  the  band  edges  where 
the  effects  of  cross  talk  are  most  severe.  With  cross  talks 
eliminated  or  ignored,  the  system  from  Sfc(n)  to  Tfc(n)  in 
Fig.  11  has  the  transfer  function  Cfc(z)  given  in  Eq.  (1). 
This  can  be  regarded  as  the  effective  channel  to  be  equal¬ 
ized  by  the  fcth  subDMT  system.  Note  that  the  cost  of 
the  filters  (Fjt(^),  Hk(z)}  is  typically  negligible  compared 
to  the  costs  of  the  DFT  matrices  involved  in  the  subDMT 
systems.  For  example  if  M  =  256  in  a  subDMT  system, 
this  costs  about  256  *  8/2  =  1,024  multiplications  at  the 
decimated  rate.  Assuming  that  all  the  DMT  systems  op¬ 
erating  at  the  decimated  rate  have  the  same  DFT  size,  the 
total  cost  of  DMT  is  1,024  multiplications  at  the  original 


rate  (i.e.,  per  sample  of  x(n)  in  Fig.  11). 

V.  CONCLUDING  REMARKS 

While  we  have  used  the  four-channel  nonuniform  system 
as  an  example  throughout,  it  is  possible  to  get  nearly  arbi¬ 
trary  tilings  of  the  frequency  axis  by  using  nonbinary  tree- 
structured  filter  banks  (wavelet  packets).  As  the  reader 
can  see,  the  purpose  of  the  paper  has  mainly  been  to  intro¬ 
duce  the  idea  rather  than  delve  into  the  details  of  many  of 
the  open  problems  that  remain  to  be  addressed.  We  know 
that  Eq.  (1)  faithfully  represents  the  transfer  function  of 
the  fcth  subchannel.  But  equally  important  are  the  cross 
talk  terms  between  the  channels,  which  have  the  form 

Ckm(z)  =  [Hk(enFm(enC(en]  ^  (2) 

where  gkm  —  gcd(njfc,nm).  Designing  the  filters  to  satisfy 
the  biorthogonality  condition  [HkFm]lgkm  =  0  for  k  ^ 
m  does  not  ensure  C/cm(z)  =  0.  Even  when  the  filters 
Fm(z )  and  Hk{z)  have  very  sharp  cutoffs,  there  is  always 
a  transition  region  of  overlap  which  creates  interchannel 
interference  and  increases  error  probability.  Finally,  even 
though  we  assign  fewer  samples  to  bands  of  lower  energy 
(e.g.,  the  band  [7T,27t]  in  Fig.  14(d)),  this  channel  also 
has  the  largest  number  of  symbols  per  second  (because  the 
decimator  and  expander  are  only  two- fold,  see  Fig.  13). 
This  will  offset,  to  some  extent,  the  advantage  obtained  by 
allocating  fewer  frequency  samples  in  that  band!  Further 
work  needs  to  be  done  to  take  into  account  these  issues, 
before  the  nonuniform  cyclic  prefix  is  deemed  a  practical 
alternative. 
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